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Abstract 

Mass Varying neutrino mechanisms were proposed to link the neutrino mass scale with dark 
energy, addressing the coincidence problem. In some scenarios this mass can present a dependence 
on the baryonic density felt by neutrinos, creating an effective neutrino mass that depends both on 
the neutrino and baryonic densities. In this article we investigate the possibility that a neutrino 
effective mass is the only flavour conversion mechanism acting in neutrino oscillation experiments. 
We present a parameterization on the environmental effects on neutrino mass that produces the 
right flavour conversion probabilities for solar and terrestrial neutrinos experiments. 
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I. INTRODUCTION 



The proposal that Mass Varying Neutrinos are coupled to dark energy [l| results in a fluid 
with negative pressure that could mimic the effects of a cosmological constant and induce a 
cosmic acceleration. The cosmological consequences of such coupling were widely addressed 



in a number of papers 

Effects on neutrino oscillation 
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were analyzed in 



different scenarios. In particular, it was proposed that the Mass Varying Neutrinos phe- 
nomenology could lead to a neutrino mass dependence with baryonic density 19| through 
non-renormalizable operators that couple the acceleron to the baryonic matter. In ref. 23 1 
some limits on the product of the effective neutrino-scalar and matter-scalar Yukawa cou- 
plings were obtained by comparing the solar neutrino and KamLAND data, assuming that 
MaVaN's mechanism plays a sub-leading role in neutrino flavour conversion. 

We present in this article a possible scenario where the neutrinos have a vanishing mass 
and mixing in vacuum. The positive oscillation indications in terrestrial and solar neutrino 
data are fully explained due to environmental effects that generate an effective mass and 



mixing 
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M 



in presence of baryonic matter. We perform an analysis of solar data 



29, 



30 



31 



35 



36 



37| , for one speciflc choice of parameters, obtaining an acceptabe solution to 



solar neutrino problem. 

In Sec. mi we introduce the general theoretical framework which we will consider in this 
paper. In Sec. IIIII we show how reactor and accelerator neutrino experiments test only the 
effective neutrino mass. In Sec. IIVI we discuss the neutrino flavour conversion probability in 
the Sun for this mechanism. In Sec. I VI I we address the problem for atmospheric neutrinos. 
Finally, in Sec. IVIIt we discuss our results and summarize our conclusions. 



II. MAVAN'S MECHANISM AND PARAMETERIZATION 



In a previous work 23|] we found limits for the product of the effective neutrino-scalar and 
matter-scalar Yukawa couplings described in [3] • But an assumption about the adiabaticity 
of the transition was made, and with this assumption we found that the new physics evoked 
always plays the role of a sub-leading effect compared to the standard oscillation scenario. 

The aim of this work is try to flnd if there is at least one combination of parameters for 
new physics that could lead to an acceptable solution to the neutrino oscillation data where 
such new physics is more then a sub- leading effect. This implies that non-adiabatic effects 
would be present in solar neutrino evolution. 

A particular case from this oscillation plus MaVaN's scenario would be the extreme 
opposite from the one investigated previously, i.e. the situation where MaVaN's physics is 
the main flavour conversion mechanism in terrestrial experiments. This is the scenario we 
investigate in this paper, where the new physics generates the neutrino mass for all terrestrial 
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indications for neutrino oscillations, and is present in all three neutrino families. 
The mass matrix in flavour basis has the form: 



M = U. 
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where Um is the mixing matrix and Mi are the mass eigenvalues related to MaVaN's effect, 
while Uy and are the mixing matrix and mass eigenvalues in vacuum. We assume that 
the mixing angles due to MaVaN's effects are constants. 

The environment effect is introduced as a dependence of the mass terms with the baryonic 
matter density with the following parameterization: 



Mi = Mo tanh A, 



P 



(2) 



3 g/cm^ 

This parameterization was chosen to reproduce two features that we want our mass matrix 
to present: 

1. a linear growth of mass with baryonic density for small values of this density. This is 
the behavior suggested in [l^, assuming a small shift in the value of A with respect 
to its ground value. 

2. a saturation of the environmental dependence of neutrino masses for large values of 
the baryonic density. 



III. REACTOR AND ACCELERATOR NEUTRINOS 

If we assume a constant Earth density in the crust, the neutrino oscillation probabilities 
can be written analytically, with the same form of the the vacuum oscillation probabilities 
in the known mass-induced oscillation scenario. The standard mass and mixing angles are 
replaced by the effective mass and mixing angles in matter. 



The positive indications for neutrino oscillation in Earth experiments |38l . l39l . l4d . |4l| tell 
us that the mixing angles that diagonalize this mass matrix in presence of Earth matter are: 

sin^ 2^23 > 0.90 (90% C.L.) 
0.4 < tan^^ia < 2.4 (95% C.L.) 

sin2 2~^i3 < 0.1 (90% C.L.) (3) 

with the following mass squared differences: 

Arhl-^ ^ 7.6 X 10"^eV^ ; Am^^ ^ 2.5 x 10"^eV^ (4) 

We want to investigate the scenario where the MaVaN's are the leading neutrino flavour 
conversion mechanism, so we assume that the vacuum mass eigenvalues are very small, 
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Am^j << 10~^ eV^, and all the flavour conversion is induced by the flrst term in right-hand 
side of Eq. [T] 

Assuming a constant Earth crust density of p ~ 3 g/cm^, we tune our parameters in eq. [2] 
in order to reproduce the above Am^. This can be achieved for instance with the following 
choice: 

Ai = ; A2 = 0.18 ; A3 = 10 

Mo = 5 X 10-2eV (5) 

With this parameterization the mass of the second family is close to a linear regime for 
the baryonic densities present at Earth, and for the crust density of p ~ 3 g/cm^ we obtain 
Arfigi = 7.9 X 10^ eV^, with a dependence. For the atmospheric neutrino scale, the third 
mass eigenvalue is already saturated for the Earth crust density, leading to Am|2 = 2.4 x 10^ 

Choosing a convenient set of mixing angles, the oscillation of terrestrial neutrino experi- 
ments are satisfactorily explained by this parameterization. 



IV. SOLAR NEUTRINOS 



For solar neutrinos, besides the environment effect in the mass matrix, the standard 
matter interaction term will have an important role in neutrino conversion. The evolution 
matrix has the form: 
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where M is the mass matrix of Eq. [TJ 

The standard interaction term would be negligible at solar neutrino production point for 
the energies of interest, indicating that the mixing matrix that diagonalizes the full evolution 
matrix would be the same in KamLAND and in the center of the Sun. 

However, the mass term decreases faster than the standard interaction one as the neutrino 
travels towards the Sun surface, and these terms become of the same order around rsun ~ 0.8 
for typical solar neutrino energies. The mixing matrix would start to feel the modiflcations 
due to the standard interaction term Vcc and the vacuum mass terms as the neutrino 
approaches the surface of the Sun. 

In Fig. [1] and Fig. [2] we can see the behavior of the eigenvalues of the evolution matrix 
and the mixing angles inside the sun, for energies of 10 MeV and 1 MeV, respectively. We 
have chosen for the vacuum parameters: 



^13 = 

Am2i 



6*23 = 6*12 = 

= 10'^ eV^ 



Am32 



2 X 10^9 eV^ 



(7) 
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In this scenario the third family would assume a constant value for its mass eigenstate 
around r ~ 0.6, while the second family would achieve the saturation value around r ~ 0.2, 
corresponding to baryonic densities of 0.4 and 30 g/cm^, respectively. 

We can estimate the survival neutrino probability assuming that the transition is adi- 
abatic until very close to the solar surface, and then the neutrino assumes the vacuum 
values of mass and mixing in an extremely non-adiabatic transition. The electron neutrino 
is created as an admixture of the two first mass eigenstates. 



>= cosOklWi > + sin I > 



where I wrote 6kl to make it explicit that the mixing angle measured in KamLAND is the 
same at Sun's interior, 6kl = (^i2{f = 0). For = 10 MeV (fig [I]) the mixing angle increases 
as the neutrino travels towards the Sun's surface due to Vcci and at r ~ 0.95 the mixing 
angle achieves the value sin ^12 = 1. Assuming this change to be adiabatic, the probability 
to have an electron neutrino at this point can be estimated to be: 

Pee = sin^ 6kl sin^ 6^12 + cos^ 6kl cos^ 6^12 = sin^ 6kl ~ 0.3 , 

and the conversion probabilities are: 

= Per ~ 0.35 

Since in an extremely non-adiabatic transition there is no conversion between flavour eigen- 
states, these are the conversion probabilities right outside the Sun. Assuming no mixing in 
vacuum, these probabilities correspond to the probabilities measured at Earth, and to the 
admixture of mass eigenstates that evolve from Sun to Earth. 

The Earth regeneration is expected to be small in this scenario. The key feature is that the 
vacuum mixing angle ^23 was made very small, while ^23 is maximum close to solar border. So 
the extremely non-adiabatic transition from the outer parts of the Sun to the vacuum leads 
to a strong production of the mass eigenstate z/3. And since we have Pee ~ -Pe^t ~ Per ~ 1/3 
the non-adiabatic transition to vacuum will lead to an equipartition of i^i, V2 and z/3, and then 
no regeneration effect at Earth. So the absence of Earth regeneration in Super-Kamiokande 
and SNO is a direct consequence of the probability Pee ~ 1/3 inside the Sun. 

Doing the same analysis for E = 1 MeV, we can estimate the probabilities the same way, 
with the difference that the mixing angle just before the non-adiabatic transition to vacuum 
is not changed by standard interaction term. Then: 

Pee = sin"^ 9kl + cos^ 9kl ~ 0.6 

and 

Pe^. = Per ~ 0.2 

Again these probabilities correspond to the probabilities at Earth if we assume non- 
adiabatic transitions and no vacuum mixing. For these energies we would expect a strong 
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FIG. 1: Mass eigenvalues and mixing angles for = 10 MeV. 



effect due to interaction with Earth matter, but since P^^ > P^.^ , the Earth effect would lead 
to a stronger conversion at night. 

We numerically calculated the solar neutrino survival probability for the values of 
tan^ 6*12 = 0.4, tan ^23 = 1 and ^13 = for the mixing angles. For the vacuum parame- 
ters, we assumed very small values for the mass eigenvalues and mixing angles. 

The result of our calculation can be seen in Fig. [31 All energy dependence present in the 
probability comes from non-adiabatic effects in neutrino evolution close to the border of the 
Sun. 

The probability obtained reproduces two main ingredients of the desired conversion to 
explain the solar neutrino data: a higher survival probability for low energy neutrinos and a 
small regeneration effect for high energy neutrinos. Besides, there is one interesting feature 
of this mechanism that is the positive day-night asymmetry for low energy solar neutrinos, 
which could be tested by Borexino. Also, this would lead to a negative winter-summer 
asymmetry in low energy solar neutrino experiments due to the difference of day and night 
duration in winter and summer. GNO recently reported a winter-summer asymmetry of 
^{W — S) = —7.6 ±8.4 SNU (~ —11%) in their full data analysis. The expected value from 
the modulation only is +2.3 SNU (+3.3%), slightly la above GNO results. 
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FIG. 2: Mass eigenvalues and mixing angles for Eiy = 1 MeV 



V. NUMERICAL ANALYSIS 



We present in this section the results of the numerical analysis of solar data for the 
parameterization discussed in the previous section. Details of our solar neutrino analysis 
have been described in previous papers 23|, |42|, |43| . We use the solar fluxes from Bahcall and 
Serenelli (2005) 4J]. In comparison with previous works, we include the new SNO data 36| . 
and included late" Borexin^ .e=nto 0- Beside., Gallex/GNO result, were split by winter 
and summer data. The solar neutrino data includes a total of 124 data points: 



1 data point for Homestake results 29| 



1 data point for SAGE results 30 1. 

2 data points for Gallex/GNO results, for winter and summer 



31 



32|. 



44 data points for Super-Kamiokande zenithal/spectral bins 331 ]. 



34 data points for SNO, phase 1 |3J]. 



38 points for SNO, phase 2 |35|. 



1 point for Borexino data 37] 



3 points for SNO, phase 3 [36 1 




Energy (MeV) 



FIG. 3: Solar neutrinos survival probability during the day (solid line) and the night (dashed line). 

We obtaine a viable solution to solar neutrino problem with the following parameters. 

tan^ Okl = 0.4 Am^^ = 7.9 x 10"^ eV^ , 

with a = 118.9. For comparison, our standard analysis provides a = 113 for the same 
parameters values. 

As mentioned before, a clear signature of this mechanism would be the day-night asym- 
metry for low energy neutrinos. For the point specified above, the Berilium line neutrinos 
would have: 

Pday = 0.62 Pnight = 0.45 , 

leading to a day-night asymmetry in Borexino oi Ajj^ = 30%. 

For the winter-summer asymmetry in low-energy experiments, this same point would 
predict: 

R summer 68.0 Rwinter 68.7, 

leading to a winter-summer asymmetry of A^s = +1%. MSW prediction for such asymme- 
try at the b.f.p. is around Aws = +4%. 



VI. ATMOSPHERIC NEUTRINOS 

In this scenario, atmospheric neutrinos would oscillate through an almost constant Am^ 
inside the Earth, with a small decrease in its value in Earth's core due to the rise of the second 
mass scale with higher densities. The possibility of environmental effects on atmospheric 
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FIG. 4: Number of oscillation lengths for atmospheric neutrinos for our scenario (continuous line), 
standard oscillation (dotted line) and previous parametrization on MaVaN's mechanism (dashed 
line). 



neutrinos were analyzed 27| with other parameterization of the matter-dependence. In that 
work the following choice: 



Amlff 



1.95 X 10" 



-0.04 



eV^ 



led to an acceptable solution to the SK atmospheric neutrino data. Here pe is the electron 



neutrino density and po = 6.02 x 10 cm~ . 

To compare our proposal with the one above, we plot in fig. H] the number of oscillation 
lengths covered by a typical atmospheric neutrino, given by: 



/A 



m^ffdL 



1 GeV 



The dashed line represents the choice in 27|, and the solid line the 6m obtained with our 
parameterization. Also presented in dotted line is the standard oscillation scenario with 
constant Am^ = 2.3 x 10"^ eV^ 

We can see that the main difference between our choice of parameterization and the 
standard scenario happens for down-going neutrinos, with intermediate energy, E^, ~ few 
GeV. There are some indications that these neutrinos are already oscillating, which would 
decrease the concordance of this model with atmospheric neutrino oscillation data. But 
due to the very similar agreement with standard oscillation in a large range of atmospheric 
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neutrino flux parameters, we believe that an overall flt would give a viable solution also to 
atmospheric neutrinos. A detailed numerical analysis would be necessary to verify this issue. 

VII. CONCLUSIONS 

We present in this article a possibility where all flavour conversion on neutrinos can 
come from environmental effects. The standard oscillation scenario is still the most elegant 
theoretical framework that explains neutrino flavour conversion, not only due to the excellent 
numerical flt to all oscillation data, but also due to the success of the research program that 
predicted, for instance, the correct signal at KamLAND based on one possible solution to the 
solar neutrino problem. The model proposed here explains such concordance a posteriori. 
However, it also make some very particular predictions on future experiments, that can 
be verifled in the close future. Borexino experiment is already taking data, and very soon 
can test a possible day-night asymmetry for low energy solar neutrinos. The author does 
not know any other model that predicts such asymmetry, making this prediction an unique 
signature of MaVaN's as the mechanism on neutrino flavour conversion. 

Acknowledgment s 

The author is thankful to R. Zukanovich-Funchal, M. C. Gonzalez-Garcia and M. Maltoni 
for fruitful discussions, and for CNPq and Fapesp for flnancial support. 



[1] R. Fardon, A. E. Nelson and N. Weiner, JCAP 0410, 005 (2004) [astro-ph/0309800l ; see also 

P. Gu, X. Wang and X. Zhang, Phys. Rev. D 68, 087301 (2003). 
[2] P. Q. Hung and H. Pas, Mod. Phys. Lett. A 20, 1209 (2005). 
[3] R. D. Peccei, Phys. Rev. D 71, 023527 (2005). 
[4] H. Li, Z. g. Dai and X. m. Zhang, Phys. Rev. D 71, 113003 (2005). 
[5] X. J. Bi, B. Feng, H. Li and X. m. Zhang, Phys. Rev. D 72, 123523 (2005). 
[6] A. W. Brookfield, C. van de Bruck, D. F. Mota and D. Tocchini-Valentini, Phys. Rev. Lett. 

96, 061301 (2006). 
[7] R. Horvat, JCAP 0601, 015 (2006). 

[8] N. Afshordi, M. Zaldarriaga and K. Kohri, Phys. Rev. D 72, 065024 (2005). 
[9] H. Li, B. Feng, J. Q. Xia and X. Zhang, Phys. Rev. D 73, 103503 (2006). 
[10] N. Weiner and K. M. Zurek, Phys. Rev. D 74, 023517 (2006). 

[11] A. W. Brookfield, C. van de Bruck, D. F. Mota and D. Tocchini-Valentini, Phys. Rev. D 73, 

083515 (2006) [Erratum-ibid. D 76, 049901 (2007)]. 
[12] R. Takahashi and M. Tanimoto, JHEP 0605, 021 (2006). 



10 



[13] A. Ringwald and L. Schrempp, JCAP 0610, 012 (2006). 

[14] G. B. Zhao, J. Q. Xia and X. M. Zhang, JCAP 0707, 010 (2007). 

[15] O. E. Bjaelde, A. W. Brookfield, C. van de Bruck, S. Hannestad, D. F. Mota, L. Schrempp 

and D. Tocchini-Valentini, JCAP 0801, 026 (2008). 
[16] K. Ichiki and Y. Y. Keum, JCAP 0806, 005 (2008). 
[17] K. Ichiki and Y. Y. Keum, JHEP 0806, 058 (2008). 
[18] O. E. Bjaelde and S. Hannestad. [aIXiv:0806.2146i [astro-ph]. 
[19] D. B. Kaplan, A. E. Nelson and N. Weiner, Phys. Rev. Lett. 93, 091801 (2004). 
[20] V. Barger, P. Huber and D. Marfatia, Phys. Rev. Lett. 95, 211802 (2005). 
[21] V. Barger, D. Marfatia and K. Whisnant, Phys. Rev. D 73, 013005 (2006). 
[22] M. Cirelh, M. C. Gonzalez-Garcia and C. Pena-Garay, Nucl. Phys. B 719, 219 (2005). 
[23] M. C. Gonzalez-Garcia, P. C. de Holanda and R. Zukanovich Funchal, Phys. Rev. D 73, 

033008 (2006). 

[24] P. H. Gu, X. J. Bi and X. m. Zhang, Eur. Phys. J. C 50, 655 (2007). 

[25] P. H. Gu, X. J. Bi, B. Feng, B. L. Young and X. Zhang, Chin. Phys. C 32, 530 (2008). 

[26] T. Schwetz and W. Winter, Phys. Lett. B 633, 557 (2006). 

[27] K. Abe et al. [Super-Kamiokande Collaboration], Phys. Rev. D 77, 052001 (2008). 

[28] K. M. Zurek, JHEP 0410, 058 (2004). 

[29] B.T. Cleveland et al., Astrophys. J. 496 (1998) 505. 

[30] R. L. Hahn, talk at XXIII International Conference on Neutrino Physics ans Astrophysics, 
New Zealand, 2008. 

[31] M. Altmann et al. [GNO COLLABORATION Collaboration], Phys. Lett. B 616, 174 (2005). 
[32] GALLEX cohaboration, Phys. Lett. B 447 (1999) 127. 

[33] Super-Kamiokande Collaboration, S. Fukuda et al., Phys. Rev. Lett. 86 (2001) 5651. 

[34] SNO Collaboration, Q.R. Ahmad et al., Phys. Rev. Lett. 87 (2001) 071301; Phys. Rev. Lett. 

89 (2002) 011301; SNO Collaboration, S.N. Ahmed et al., Phys. Rev. Lett. 92 (2004), 181301. 
[35] B. Aharmim et al. [SNO Collaboration], Phys. Rev. C 72, 055502 (2005). 
[36] B. Aharmim et al. [SNO Collaboration], Phys. Rev. Lett. 101, 111301 (2008). 
[37] C. Arpesella et al. [The Borexino Collaboration], Phys. Rev. Lett. 101, 091302 (2008). 
[38] T. Araki et al. [KamLAND Collaboration], Phys. Rev. Lett. 94 (2005) 081801; S. Abe et al. 

[KamLAND Collaboration], Phys. Rev. Lett. 100 (2008) 221803; KamLAND data and 

results are available at 

http:/ /www.awa.tohoku. ac.jp/KamLAND / chi2map_3rdresult / chi2map.htniil 
[39] M. Apollonio et al. [CHOOZ Collaboration], Phys. Lett. B 466 (1999) 415. M. Apollonio et 

al. [CHOOZ Collaboration], Eur. Phys. J. C 27 (2003) 331. 
[40] M. H. Ahn et al. [K2K Collaboration], Phys. Rev. D 74 (2006) 072003. 

[41] D. G. Michael et al. [MINOS Collaboration], Phys. Rev. Lett. 97 (2006) 191801. P. Adamson 
et al. [MINOS Collaboration], Phys. Rev. D 77, 072002 (2008). H. Gallagher, talk at XXIII 



11 



International Conference on Neutrino Physics ans Astrophysics, New Zealand, 2008. 
[42] J.N. Bahcall, M.C. Gonzalez-Garcia and C. Pena-Garay, JHEP 02 (2003) 009 

|hep-ph/0212147] ; J.N. Bahcall and C. Pefia-Garay, JHEP. 11 (2003) 004. 
[43] P. C. de Holanda and A. Y. Smirnov, Astropart. Phys. 21, 287 (2004). 
[44] J. N. Bahcall and A. M. Serenehi, Astrophys. J. 626, 530 (2005). 



12 



